Non-ergodic dynamics of the extended anisotropic Heisenberg chain 



Evgeny PlekhanovQ Adolfo Avella0 and Ferdinando Mancin^ 
Dipartimento di Fisica "E.R. Caianiello" - Unita CNISM di Salerno 
Universitd degli Studi di Salerno, I-84O8I Baronissi (SA), Italy 

The issue of ergodicity is often underestimated. The presence of zero- frequency excitations in 
bosonic Green's functions determine the appearance of zero-frequency momentum-dependent quan- 
tities in correlation functions. The impHcit dependence of matrix elements make such quantities 
also relevant in the computation of susceptibilities. Consequently, the correct determination of these 
quantities is of great relevance and the well-established practice of fixing them by assuming the er- 
godicity of the dynamics is quite questionable as it is not justifiable a priori by no means. In this 
manuscript, we have investigated the ergodicity of the dynamics of the z-component of the spin 
in the ID Heisenberg model with anisotropic nearest-neighbor and isotropic next-nearest-neighbor 
interactions. We have obtained the zero-temperature phase diagram in the thermodynamic limit by 
extrapolating Exact and Lanczos diagonalization results computed on chains with sizes L = 6 26. 
Two distinct non-ergodic regions have been found: one for J'/Jz < 0.3 and \J±\/Jz < 1 and another 
for J'/Jz ^ 0.25 and \ J±_\/ J^ — 1. On the contrary, finite-size scaling of T 7^ results, obtained by 
means of Exact diagonalization on chains with sizes L = 4 18, indicates an ergodic behavior of 
dynamics in the whole range of parameters. 



I. INTRODUCTION 

In order to study a physical system (i.e., to analyze its 
properties and its responses to external perturbations), 
we have to compute correlation functions C and suscep- 
tibilities X (retarded Green's functions) of any operator 
Tp relevant to the dynamics under investigation 
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According to the well-known relations existing between 
casual (C), retarded {R) Green's functions and correla- 
tion functions 



where (• ■ ■ ) stands for the average in some statistical en- 
semble, i = (i, t) for both spatial i and temporal t co- 
ordinates. The operator ip is taken in the Heisenberg 
picture. 

In principle, one can compute the complete set of eigen- 
states |n) and eigenvalues En of the Hamiltonian describ- 
ing the system under study. This knowledge makes pos- 
sible to compute the Fourier transform in frequency of 
any correlation function as follows 

C(i,j,c^) = 2^<5(u;)r(i,j)+ 

Y E e-'3^"(n|V(i)|m)(m|V/(j)|n),5[c^ + (i?„-i?™)] 

where the zero-frequency function r(i, j) is defined 
r(iJ) = | E e-'3^"(nl^(i)|m)(m|^t(j)|„). (1.3) 

r explicitly appears in the expression of the correlation 
function when the field tp is boson-like, that is, ip is con- 
stituted of an even number of original electronic (i.e., 
fermionic) operatorsiiSi. Hereafter, we focus on such a 
case. 
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3[G^(k,c.)] = tanh ( ^ ) 5[G^'(k,^)] 
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the zero-frequency excitations do not contribute explic- 
itly to the imaginary part of the retarded Green's func- 
tions and, consequently, F does not explicitly appear in 
the expressions of susceptibilities. At any rate, suscepti- 
bilities retain an implicit dependence on F through the 
matrix elements (n|';/'(i)|m). Then, the right procedure 
to compute both correlation functions and susceptibili- 
ties is clearly the one that starts from the causal Green's 
function, which is the only Green's function that explic- 
itly depends on F. 

It is worth noticing that the value of F dramati- 
cally affects the values of directly measurable quantities 
(e.g., compressibility, specific heat, magnetic susceptibil- 
ity, . . . ) through the values of correlation functions and 
susceptibilities. According to this, whenever it is possi- 
ble, F should be exactly calculated case by case. 

If we do not have access to the complete set of eigen- 
states and eigenvalues of the system, which is the rule 
in the most interesting cases, we have to compute cor- 
relation functions and susceptibilities within some, often 
approximated, analytical framework. Now, since no ana- 
lytical tool can easily determine F (e.g., the equations of 
motion cannot be used to fix F as it is constant in time) , 
one usually assumes the ergodicity of the dynamics of tp 
and simply substitutes F by its ergodic value: 



F-niJ)-(V'(i))(V'Mj)) 



(1.7) 



Unfortunately, this procedure cannot be justified a priori 
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(i.e., without computing F through its definition H1.3|) ') by 
absolutely no means. The existence of just one integral of 
motion divides the phase space into separate subspaces 
not connected by the dynamics. This latter, in turn, 
becomes non ergodic: time averages give different results 
with respect to ensemble averages. 

The lack of ergodicity has sizable effects. One of its 
manifestations is the well-known difference between the 
static isolated (or Kubo) susceptibility and the isother- 
mal one^. The dynamics of an operator ip, in a system 
described by the Hamiltonian H, is said to be ergodic iS 



hm - [ (^Vt(t))dt-(V')(^^), 



(1.8) 



where (• • • ) stands, for instance, for the canonical average 
at temperature T 



(1.9) 



with /3 = y and Z = Tr(e"''^). That is, the dynamics 
of ip is ergodic if, during its time evolution, it has non- 
zero matrix elements only between states within a zero- 
volume region of the phase space of the system. Only 
in this case, the equilibrium averages (l.h.s. of (|1.8|) ') are 
equal to the ensemble averages (r.h.s. of (|1.8|l ). which 
are much easier to compute. Now, one can show^ that 
the difference between the static (or Kubo) susceptibility 



a vanishing ■0-7i coupling in the Hamiltonian. For in- 
stance, the reader can imagine a ferromagnetic ground 
state with finite spontaneous magnetization driven by a 
vanishing external magnetic field. 

Many systems exhibit non-ergodic dynamics. Over the 
last years, we have been studyingi*^ some of them ana- 
lytically by means of the Composite Operator Method^-—: 
two-site Hubbard model, three-site Heisenberg model, 
tight-binding model in magnetic field, double-exchange 
model, extended Hubbard model in the ionic limit. Re- 
cently, we have started an investigation of such models 
also by means of numerical techniques: Lanczos and Ex- 
act Diagonalization. Along this latter line, in the present 
paper, we calculate F for the z-component of the spin at 
site i, S? , in the one-dimensional anisotropic extended 
Heisenberg model (see next section) and show that it 
takes non-ergodic values in some regions of the parame- 
ter space of the model both for finite sizes and in the bulk 
limit. This analysis is the natural continuation of that 
reported in Ref. |a, where we have analyzed the standard 
(isotropic and non-extended) Heisenberg model at zero 
temperature and found a non-ergodic behavior of F. 



II. DEFINITIONS AND METHOD 

We have studied the ergodicity of dynamics of the oper- 
ator S? in the ID anisotropic extended Heisenberg model 
described by the following Hamiltonian: 



X(0)= lini^[-i0(i)([V',Vt(i)])] 

= / {tpip^iX}) d\ - 13 lim - 
Jo 



^oo T 



(1.10) 



and the isothermal susceptibility 



^ SH 



(V'V^(iA))dA-/3(^)(V/), (1.11) 



is just proportional to the difference between the l.h.s. 
and the r.h.s. of (|1.8|l 



- X(0) = /3 



lim — 



(1.12) 



In the above, J- stands for the Fourier transform and Ti. 
for the external field coupled to ip (e.g., an external mag- 
netic field in the case of magnetization) . The second line 
of (|1.10(l comes from the relation between susceptibility 
and relaxation function. F (as defined in (jl.3|l l is exactly 
equal to the l.h.s. of (|1.8I) . 

Thus, it is clear that issue of ergodicity is relevant to 
linear response too. It is also worth noting that non- 
ergodicity at zero temperature is just the result of a de- 
generacy in the ground state, which would be lifted by 



i 

+ ^^(^r^r+i + sfsf^,) + SiSi+2, (2.1) 



where S^, Sf and S? are the x, y and z components of the 
spin-1/2 at site i, respectively. The model (|2.1|l is taken 
on a linear chain with periodic boundary conditions. We 
take the interaction term parameterized with Jz ferro- 
magnetic { Jz > 0) and the next-nearest-neighbor inter- 
action term, which is parameterized with J', isotropic. 
In order to frustrate ferromagnetism, we have considered 
only the case with J' > 0, that is, with an antiferromag- 
nctic coupling between next-nearest neighbors. Accord- 
ing to this, only chains with even number of sites have 
been studied in order to avoid topological frustration that 
would be absent in the thermodynamic limit. Since it is 
possible to exactly map all results obtained for Jj^ > to 
those for Jj_ < by means of a simple canonical transfor- 
mation, we have limited our study only to positive values 

of Jj^. 

We have found some numerical2*2i2ii2iii and few an- 
alytical studiesiSii^ of the model 12.111 in the literature, 
mainly in the antiferromagnetic state (J^ < 0). Only one 
fixed point {J± = Jz and J' = Jz/2) is known to pos- 
sess an analytic ground statai^. In the rest of the phase 
diagram, the model (|2.1|l seems not to be integrabl o^^'^^ . 



We have numerically diagonalized the Hamiltonian 
(|2.1|) for chains of size L ranging between 6 and 18 by 
means of Exact Diagonalization (ED) (divide and con- 
quer algorithm) and for chains of size L ranging between 
20 and 26 by means of Lanczos Diagonalization (LD). We 
have systematically taken into account translational sym- 
metry and classified the eigenstates by the average value 
of S*^ = Tlii^i^ which is a conserved quantity. When- 
ever we have used ED, all eigenvalues and eigenvectors 
of (|2.1|) have been calculated up to machine precision 
and, therefore, we have been able to determine the ex- 
act dynamics of the system for all temperatures. On the 
contrary, when we have used LD, we have been limited 
to the zero-temperature case since only the ground state 
can be considered exact in LD. 

As already discussed above, the dynamics of an oper- 
ator (e.g., 5?) is ergodic whenever (|1.8() is satisfied, or 
equivalently, H1.3|) is equal to its ergodic value: 



n.ni 



(2.2) 

The dynamics of a finite system is hardly ergodic, since 
and (|2.2|) unlikely coincide. In the thermodynamic 
limit, the sums in H1.3|l and H2.2() become series and no 
conclusion can be drawn a priori. Since we have diago- 
nalized the Hamiltonian (|2.1|) numerically (i.e., only for 
finite systems) and since L — > oo is the most interest- 
ing case, we have analyzed our results through finite-size 
scaling in order to speculate on the properties of the bulk 
system. 



III. ZERO-TEMPERATURE RESULTS 

If the ground state of H2.1|) is A-fold degenerate then, 
at r = 0, (|1.3|1 and (|2.2|1 read as follows: 




perg 



respectively. 

Thanks to the translational invariance enjoined by the 
system {S?) is independent of i and proportional to the 
z-component of the total spin operator average {S^^^). It 
is easy to show that, even if there is a finite magnetic mo- 
ment per site in any of such N degenerate ground states, 
(Sf ) at T = is always zero in absence of magnetic field. 
Indeed, if a ground state with non-zero (5'f ) = M exists, 
also another ground state with {S() — —M exists. Thus, 
at zero temperature, F'^''^ is always zero and the only 
quantity of interest is F. A finite value of this latter im- 
plies non-ergodicity. Obviously, if = 1 then both val- 
ues coincide. Therefore, a non-ergodic phase corresponds 
to degenerate ground states with finite magnetization. 
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FIG. 1; Zero-temperature ergodicity phase diagram in the 
J' — Jx plane. Due to the symmetry of the Hamiltonian 
only the upper half is shown (see in the text). Only two 
ergodic phases (E-I and E-II) have been found in the reported 
parameter space. The others are either non-ergodic (NE-I and 
NE-II) or impossible to conclusively analyze (W). The latter 
phase might shrink to a transition line in the bulk limit. 



In the studied range of coupling constants (see Fig. ^ 
we have found two non-ergodic phases (NE-I and NE-II), 
two ergodic ones (E-I and E-II) and a weird phase (W). 
Our computational facilities limit the range of chain sizes 
that we can analyze such that we could not establish, by 
means of finite-size scaling, whether the weird phase (W) 
is ergodic or not. In the non-ergodic phases (NE-I and 
NE-II), we were able not only to perform the finite-size 
scaling, but also to write down an analytic expression for 
F as a function of the chain size L. The weird phase 
(W) has exhibited a strong dependence of the ground 
state upon the particular values of the couplings. On 
the contrary, the other phases exhibit ground states that 
are independent of the particular values of the coupling 
constants. 

Ergodicity of the standard Heisenberg model {J' — 
and J\^ — Jz)aXT = Q has already been studied in Ref.0: 
the dynamics is non-ergodic for ferromagnetic coupling 
( J_i = Jz < 0) as the system has a L -I- 1 degenerate 
ground state 



F 
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12 



1 



(3.1) 



It is clear from (|3.1|l that F remains non-ergodic also 
in the thermodynamic limit. This point (J' = and 
Jj_ = Jz) becomes a line in our phase diagram and is 
denoted as NE-I (see Fig. In fact, the next-nearest- 
neighbor interaction J' may frustrate (J' > 0) or favor 
(J' < 0) the ferromagnetism. In the latter case, the 
ground state remains unchanged for any value of J' < 0. 
Therefore, we expect the line denoting the phase NE-I 
to extend also to negative J'. If, on the contrary, J' 
is positive and large enough to frustrate the system in 
such a way that the ground state loses its ferromagnetic 
character, the ergodicity is restored. This occurs at a 
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FIG. 2: Finite-size scaling in the case of T = for different points in the phase diagram of Fig. Symbols on panel a): 
+ corresponds to (NE-II), ▲ corresponds to (NE-I), ■ corresponds to (E-I) and Q to (E-II) regions of Fig. respectively. 
On panel b) different examples from (W) region are shown. Hamiltonian couplings are shown in the legend. All energies are 
expressed in units of Jz. 



finite critical J' ^ 0.25 Jz- For values of J' larger than 
the critical one, we find a non-degenerate ground state 
with (Sf,,) = 0. 

If J± 7^ Jz the rotational invariance is broken so that 
states with the same (Sj^t), but different (S^^^), are not 
degenerate anymore. In the non-ergodic region (NE-II) 
of the phase diagram (see Fig.P), the ground state is just 
doubly degenerate (not L + 1 degenerate as in (NE-I)): 
one ground state corresponds to a configuration with all 
spins up and the other to a configuration with all spins 
down. Hence, the value of T in this phase is 1/4 and 
does not depend neither on the Hamiltonian couplings 
nor on the number of sites in the chain. It is clear that 
also this phase extends to negative values of J'. This 
kind of ground state stands the frustration introduced 
by next-nearest-neighbor interaction up to J' ~ O.SJ^ 
(see Fig.[Tl). 

The ergodic region (E-I) of the phase diagram (see 
Fig. ^ has r = for all sizes of the system and val- 
ues of the couplings: the unique ground state belongs to 
the sector with (S'fot) — 0. On the contrary, the other er- 
godic phase (E-II) has non-zero values of F for values of L 
not multiples of four. The ground state in this phase has 
average total spin equal to one and, therefore, F = 
We obviously conclude that (E-II) phase is ergodic in the 
thermodynamic limit. 

The values of F in these four phases (NE-I, NE-II, E-I 
and E-II) exhibit perfect finite-size scaling as shown in 
Fig. Ek) • This has allowed us to make definite statements 
also in the thermodynamic limit. 

The weird phase ( W) (see Fig. ^ is characterized by a 
quite strong size dependence, as shown in Fig.j^b) where 
a tentative finite-size scaling of F in the different points 
of the phase is presented. This region manifests a di- 
verging finite-size scaling within the range of sizes we 
were able to handle. In this case, the behavior of F as a 
function of L strongly depends on the particular choice of 



the Hamiltonian couplings and is highly non monotonous 
when increasing L, according to the strong dependence 
on L of (Sfot) in the ground state. In this critical re- 
gion the eigenvalues of (|2.1I) present many level crossings, 
which means that the maximum value of L we were able 
to reach {Lmax = 26) is not large enough to perform a 
sensible finite-size scaling analysis. However, we expect 
that this phase becomes ergodic in the thermodynamic 
limit, although still different from the ergodic phases E-I 
and E-II. 

We can summarize our findings in the thermodynamic 
limit at zero temperature as follows: 

' if = ±Jz and J' < 0.25 

if I J± I and J' < 0.3 J, ^^^^ 
111 in the weird phase (W) (see Fig. 
otherwise 

V 



IV. FINITE-TEMPERATURE RESULTS 

At zero temperature, the degeneracy of the ground 
state is clearly the source of non-ergodicity. At finite 
temperatures, all states contribute to the dynamics of 
the system and we have to search for some other source 
of non-ergodicity, if any. In fact, when the sum in \1.'6\ 
becomes an integral, the degree of degeneracy of a state 
with energy e is given by the value of the density of states 
at £. Then, the result of the integration is hard to guess. 
Once again, we have used finite-size scaling in order to 
reveal the behavior in the bulk system. Ergodic value of 
F should now be calculated through the thermal average: 



re.,^ _^^-/5iJ„^^^|^.|^^ • (4-1) 
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FIG. 3: Finite-size scaling of the finite-temperature results for different representative points on tlie pfiase diagram of Fig.0 a) 
NE-I pliase. b) NE-II pliase. c), d) two sample points from the weird (W) phase, which at T = show quite distinct finite-size 
scaling and at T > behave very similarly. All energies are measured in units of Jz- AH the data suggest that the bulk system 
at T > is always ergodic. 



It follows that r'^'"^ is just zero for all Hamiltonian cou- 
plings. 

We have computed finite-size scaling at several finite 
temperatures using ED. In Fig. |3| four representative 
points from the phase diagram in Fig. ^ are shown: the 
plots from panels a) and b) belong to the former NE- 
I and NE-II phases, while those from panels c) and d) 
belong to the W one. Finite-size scaling for the ergodic 
phases (E-I and E-II) gives F identically zero and is not 
reported. Already at low temperatures (T ~ 0.08 Jz) and 
small sizes L < 18, our calculations show a very clear 
indication that, in the thermodynamic limit, F is zero 
in the whole phase diagram. In addition, we have found 
that F approaches zero nearly proportional to L^^ for 
finite systems. 

These results can be analyzed from another point of 
view. As it is stated by Suzuki^, necessary and sufficient 
condition for the ergodicity of the dynamics of an oper- 
ator A under the Hamiltonian H is the orthogonality of 



A to all the integrals of motion of the system 

{A,W)^Q, (4.2) 

where {i?*} represents the complete set of integrals of 
motion of the Hamiltonian H. {H'^} form a linear space 
with metrics given by the normalized density matrix 

{H\ W) = 1 ^ e-P^^Hl,Hl, ^ {WW), (4.3) 
fe 

where iJ* and are two integrals of motion and the 
sum is extended on all microscopic states of the system. 
Then, for a given operator A, to be orthogonal to all in- 
tegrals of motion of H means either (i) to be off-diagonal 
with respect to H or (ii) to have a negligible number 
of diagonal matrix elements so that the scalar products 
(|4.2|) would go to zero in the thermodynamic limit. It fol- 
lows that the dynamics of an integral of motion is surely 
non-ergodic, since it is surely not orthogonal to itself. 

For finite systems, the response is practically always 
non-ergodic as it is unlikely that H1.3|) is exactly zero. In 
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fact, for any fixed size, F is non zero (see Fig.O. Naively, 
one could expect that S? has non-ergodic dynamics also 
at finite temperatures since its average is proportional to 
the one of S^^^, which is an integral of motion. However, 
our numerical studies suggest that case (ii) applies to this 
operator. It is worth noticing that this type of behavior 
can be strongly dependent on the dimensionality of the 
system. For instance, we can expect a non-ergodic dy- 
namics also at finite temperatures in higher dimensions, 
where a finite critical temperature for the transition to 
the ferromagnetic phase exists. 

As it results from our analysis, the question of ergodic- 
ity of a spin system is closely connected to the presence of 
a finite magnetic moment per site. It might seem, there- 
fore, that since we are dealing with a one-dimensional 
system, we could apply the Mermin- Wagner theorem^^ 
and exclude the possibility to have finite magnetization 
at T ^ 0. Actually, our model is not isotropic so one 
of the conditions of Mermin- Wagner theorem is not met. 
At any rate, even along the line J± — Jz the proof of the 
above theorem involves the introduction of an infinitesi- 
mally small magnetic field that would split the degener- 
acy and make ergodic the system dynamics. 

V. CONCLUSIONS 

Summarizing, we have studied the ID anisotropic 
Heisenberg model with next-nearest-neighbor exchange 



interaction and the ergodicity of the operator S'f by 
means of the Lanczos and Exact Diagonalization tech- 
niques. At zero temperature, we have constructed the 
ergodicity phase diagram and found two distinct non- 
crgodic regions in the thermodynamic limit. The borders 
of these regions and the corresponding non-ergodic values 
of F have been determined by means of a finite-size scal- 
ing analysis of the Lanczos data. The results of this anal- 
ysis (i.e., critical couphngs and non-ergodic F values, see 
(j^l) can be used by analytic methods that aim at com- 
puting physical properties and response functions of this 
system. Our results for finite temperatures suggest the 
ergodicity of the dynamics in the thermodynamic limit in 
the whole range of couplings. It implies, on the basis of 
Suzuki theoreni2,^ that the number of diagonal matrix el- 
ements of 5? in the energy representation is negligible in 
the thermodynamic limit. A finite-size system, instead, 
remains non ergodic and the finite-L corrections to the 
bulk value of F are proportional to the inverse volume of 
the system (i.e., to L^^)- 
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